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Abstract. We propose a new method to obtain approximate solutions for the 
Schrodinger equation with an arbitrary potential that possesses bound states. This 
method, relying on the auxiliary field technique, allows in many cases to find analytical 
solutions. It offers a convenient way to study the qualitative features of the energy 
spectrum of bound states in any potential. In particular, we illustrate our method by 
solving the case of central potentials with power-law form and with logarithmic form. 
For these types of potentials, we propose very accurate analytical energy formulae 
which improve a lot the corresponding formulae that can be found in literature. 
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1. Introduction 

From the early years of quantum mechanics, searching for analytical solutions of the 
Schrddinger equation has always been a subject of intense investigations, in particular 
when the Hamiltonian possesses bound states. This is a challenging mathematical 
problem often involving complex differential equations, and several methods have 
been proposed to find approximate analytical solutions: WKB method, semiclassical 
treatment, variational method, perturbation theory, etc. More details can be found 
in [1] for example, in which a wide range of analytical results concerning bound states 
are presented. Nowadays, numerical calculations allow us to solve eigenequations very 
accurately. But, finding analytical solutions is still of great interest since it makes 
explicitly appear the dependence of the eigenenergies and wave functions on the various 
parameters of the model. Such a property is particularly useful when one tries to fit 
the parameters of a model to experimental data for example. In the present work, we 
propose a new way to obtain analytical approximate solutions of Schrodinger's equations 
for bound state problems. Our method is based on the auxiliary field technique. 

Auxiliary fields, also known as einbein fields, appear in various domains of 
theoretical physics. They are commonly used to get rid of the square roots appearing 
in relativistic systems; let us mention the Nambu-Goto Lagrangian which is replaced by 
the Polyakov one in string theory for example [2] . They also appear in supersymmetric 
field theories [3] and in hadronic physics [I]. Let us briefly present the auxiliary fields 
with the simple example of a free relativistic particle, that is the Lagrangian 

C = —m Vr, (1) 

where x is the world-velocity. An auxiliary field \i can be introduced in this expression 
to get rid of the square root 

x x 2 m 2 

This last expression is formally simpler than Lagrangian (CQ), but it is equivalent if 
the equations of motion of the auxiliary field are considered. The auxiliary field is 
determined by the condition of extremum 

^)U =0^ = -^, (3) 

and if we substitute it into the corresponding Lagrangian (T5]), we recover the original 
one ([1]): C(fio) = C This is a general feature of auxiliary fields: Their equation of 
motion are not dynamical but rather algebraic relations that express them in terms of 
the other variables of the problem. 

We see that the introduction of auxiliary fields can notably simplify the formalism, 
and thus the computations. It was pointed out in [5] that relativistic Hamiltonians can 
be rewritten as apparently nonrelativistic ones by such a method [compare equations ([I]) 
and (T2]) for example]. We previously showed that auxiliary fields allow to express a 
linear potential, which is widely used in hadronic physics, as an expression involving an 
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auxiliary field and a harmonic oscillator By doing this, analytical mass spectra can 
be computed assuming that the auxiliary field is not an operator but a real number. That 
is why approximate results are obtained. The question is: Could any arbitrary potential 
be expressed in terms of an equivalent expression involving an auxiliary field and another 
potential for which an analytical solution is known? As we will show in section [2} the 
answer is yes, and a general algorithm for analytically solving the Schrddinger equation 
with any potential can be proposed. In sections [3] and HI we illustrate our general 
method by computing analytical energy spectra for central potentials with power-law 
form and with logarithmic form. As our formulae are approximate, it is of interest 
to try to maximally improve them. This is done in section El Finally, we draw some 
conclusions in section [7J 



2. General method 

We begin with a standard Hamiltonian containing a kinetic energy term T(p 2 ) and a 
local central potential V(r) 

H = T(p 2 ) + V(r), (4) 

whose eigenvalues and eigenstates are a priori not analytically known. On the contrary, 
let us assume that the Hamiltonian with potential P(r) replacing V(r), namely 
Ha(p) = T(p 2 ) + pP(r), where p is a real parameter, possesses bound states with 
an analytical spectrum. Thus, we explicitly know the solution of the corresponding 
Schrddinger equation 

H A {p) = E A {p) . (5) 

Starting from the eigenequation (jHJ), we propose hereafter a way to compute an analytical 
approximate solution of the eigenequation 

H\V) = E\V), (6) 

H being the arbitrary Hamiltonian (jlj). The procedure consists in three steps: 

(i) Build the function K(r) and the operator p by the definitions 

so that the operator r is given by the inverse function 

r = K-\p). (8) 

(ii) Introduce the Hamiltonian H(p) through the definition 

H(p) = T(p 2 ) + pP(r) + V [K-\p)] -pP [K-\p)] . (9) 

p is known as an auxiliary field. Because of fl5]), the spectrum of H(p) is analytically 
known if the auxiliary field is assumed to be a number, and not an operator. The 
eigenstates are |\& (/?)}, while the eigenenergies E(p) are obtained from E A (p) simply 
by adding a p-dependent constant term 

E{p) = E A (p) + V [K-\p)] - pP [K-\p)] . (10) 
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(iii) The physical spectrum is finally given by the eigenstates \^f(p )) and the 
eigenenergies E(p ), with p minimizing the total energy ffTUj) . i.e. 
dE{p) 



dp 



= 0. (11) 

P=P0 

It is now time to justify the procedure given above. The Hamiltonian H is replaced 
by H(p), given by Actually, one replaces in H the potential V(r), for which the 
energies are not analytical, by the potential P(r) for which they are, but via an effective 
potential V(r, p) = pP{r) + V [K~ l (p)] — p P [K~ l (p)] depending on the auxiliary field. 
Such an auxiliary field is nothing else but a mathematical trick to simplify the formalism. 
In the present case for example, V(r) is "simulated" by a function of p and of an 
analytically solvable potential P[r). The key point is that the auxiliary field, when 
properly eliminated as an operator, should lead to Hamiltonian (T4j). It can be checked 
that 

5V(r,p) = p = K(r). (12) 
p=p 

Moreover, V(r, K(r)) = V(r). This means that Hamiltonians (jl]) and (Q are actually 
equivalent up to the elimination of the auxiliary field as an operator, the auxiliary field 
being then defined by (J7J). 

The approximation in our method is the following: If p is considered as a number 
in Hamiltonian ([9]), the corresponding eigenequation can be analytically solved, since 
only P(r) is involved. Then, an energy spectrum of the form fflOl) . which depends on the 
auxiliary field, is obtained. Obviously, the physical value of p for a given eigenenergy is 
the one which ensures this energy to be minimal. This is expressed in ( TTTj) . 

The Hellmann-Feynman theorem [Tj states that 

«M = <,(,,, «M ,*(,)). (13) 

Using this relation and defining the quantity tq by po = K(ro), one can show that 

(*(po)|P(r)|*(pd)> = P(r ). (14) 

This means that r is a kind of "average point" for the potential P(r). If functions P(r) 
and V(r) are smooth enough, one could expect that the average point for the function 
K{r) does not differ significantly from r$. So we can assume that 

(*(po)|p|^(po))«^(r )=po. (15) 

Thus, we can expect that the optimal value po should be close to the average value 
(\P (po) I p |^(po))) with p given by <^ [8J. Our method could actually be considered 
as a "mean field approximation" with respect to a particular auxiliary field which is 
introduced to simplify the calculations. This will be checked in the following. 

It is possible to obtain an analytical estimation of the accuracy of the method. 
With the notations introduced above, let us compute (ty(p )\H(p ) — H\^(p )). Using 
®, (ED and flUD, we find 

E(p ) - (*(p )!if|*(po)> = V(r„) - (*(pd)|V(r)|*(pD)>. (16) 
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The r.h.s. of this equation is the difference between the value of potential V computed 
at the "average point" r and the average of this potential for the trial state |\l/(po))- 
In some favorable cases |9] (the trial state is a ground state for instance), E < 
(p )\H\^/ (p )) , where E is the exact eigenvalue of H. Since this inequality cannot be 
always guaranteed, we just can write in general 



In principle, all the quantities present in this formula, except E obviously, can be 
analytically computed. So an estimation of the error on the absolute eigenvalue can 
be obtained. Depending on the sign of the r.h.s. of (fTTI) . E(p ) can be an upper or 
a lower bound of E. An estimation of the relative error is obtained by computing 
(V(r ) - (^(po)\V(r)\^(po)))/\E(p )\. Numerical evaluations of formula (JTTj) will be 
given in the following. 

Formally, we thus have proposed an algorithm which allows to compute approximate 
analytical eigenvalues and eigenstates of an arbitrary Hamiltonian containing the 
potential V(r), provided that an analytical solution of the Hamiltonian is known for 
a particular potential P(r). This method leads to approximate solutions because of 
the presence of the auxiliary field, but, as we will see, it is a very convenient tool to 
understand the qualitative features of various energy spectra. 

Before giving explicit examples, three points should be stressed: 

• One easily checks that the method is exact for V(r) = P(r), with p — 1 trivially. 
One thus expects that, the closer is P(r) to V(r), the better the approximation is. 

• Generally, the procedure is not variational, so that we have no certainties concerning 
the position of the approximate energies with respect to the exact ones. 

• The method is not free of technical difficulties. Although the presented algorithm 
is completely general and in principle applicable to any situation, obtaining an 
analytical solution is not always possible or can lead to cumbersome calculations. 
In fact the crucial difficulty concerns the inversion of the function K(r) [formula (jSJ)] 
and/or the determination of the best auxiliary field po [formula (jlip ]. However, 
these difficulties can be overcome in many relevant cases and the corresponding 
equations can be fortunately solved. 

We can also notice that the method which is proposed here remains valid not only 
for central potentials, but also for general quantum mechanical systems (one- and two- 
dimensional ones for example). 

3. The case P(r) = r 2 
3.1. Power-law potentials 

Let us first illustrate the general procedure presented in the previous section by solving 
the eigenequation associated to the nonrelativistic Hamiltonian 




(17) 



H 



2m 



+ sgn(A) ar , 



(18) 
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with A 7^ 0, sgn(A) = |A|/A, and a > in order that this Hamiltonian possesses bound 
states. 

Two power-law potentials are standard in quantum mechanics formalism: The 
harmonic and Coulomb potentials because analytical solutions are known whatever the 
radial and orbital quantum numbers. This should thus be of interest to take P(r) = r 2 
or — 1/r in our auxiliary field method. In this section, we will consider the case of 
the harmonic oscillator. We recall that the spectrum of the Schrddinger equation with 
potential 

V H {r) = ^mu 2 r 2 (19) 

is given by (see for example [1] problem 66]) 

R H (r) = N n £ / e - mwr2/2 L £ n +1/2 (mujr 2 ), E H = oo(2n + £ + 3/2), (20) 

where N n e is a normalization coefficient and R H (r) the radial wave function (n = 
0, 1, 2, . . . is the radial quantum number and i = 0, 1, 2, . . . is the orbital quantum 
number). 

Taking P(r) = r 2 , ([7]) leads to 

P=^-^r=K-^ (2D 

We notice that for A ^ 2, K (p) is always well-defined since a|A| > by definition. 
The case A = 2 is actually of no interest since in this case p = 1 trivially, and the results 
are exactly known. Then, following (Q, one obtains after some algebra 

H(p) = ^ + Pr 2 + 2 -^(2p)^(a\X\)^. (22) 

We can easily check that Hamiltonian (122]) is fully equivalent to Hamiltonian f|T8|) if 
p = p as given by ( 12TI) . In order to go further in the computation, we now consider the 
auxiliary field p as a real parameter. As we already mentioned, this is at this place that 
the approximation is introduced. This procedure allows to obtain an analytical energy 
formula which reads 

E H (p) = \[^(2n + l + 3/2) + £zi(2p)xb( |A|)^. (23) 
V m 2A 

The wave function is also readily known from (1201) in which uj is replaced by ^2p/m. 
We will not write explicitly the wave functions in the following, since they can readily 
be obtained once the energy is known. 

The last step is the minimization of the energy with respect to the auxiliary field 
p. One finds 

2(A-2) 

"2n + £ + 3/2 



1 

Po 



2 



id 



4 



A+2 HA|)^, (24) 



and then the final form for the energy spectrum is 

2A 



£"M=(^^^^(2« + ^3/2)^. (25) 



m x + 2 
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It is interesting to notice that 

E H (po)\ x=2 = \/^(2n + £ + 3/2), (26) 




EH (Po)\x=-i ~ ~ 2 {2n + £ + 3/2f (27) 

The formula is exact for A = 2 as expected, since in this case V(r) = P(r). If, in 
the case A = —1, we substitute 2n + 3/2 by n + 1, we recover the exact result for the 
Coulomb case. Of course, the presence of a term 2n + I + 3/2 is a "remembrance" of 
the original function P(r) that we employed to treat the auxiliary field. The fact that a 
potential with a very strong increase for large r can mimic the asymptotic behavior of 
the energies for a potential that tends to zero for large r is astonishing and shows the 
power of our method to get at least the good asymptotic properties. We remark also 
that formula ff25l) is not defined when A = —2. Moreover, the energy becomes positive 
when A < —2. We recover actually a well-known result: In an attractive potential with 
a singularity like l/r n with n > 2, the existence of bound states is not automatically 
guaranteed [10J. 

We mentioned in the previous section that po was expected to be close to the average 
value of p. The average should here be computed with the wave functions R H (r) given 
by ( l20l) in which oj = ^j2~pofrn. This point can be explicitly checked in the particular 
case we are dealing with. Indeed, relation (1211) shows that 

W = ^l( r «) fs ^^) i f 2 . ( 28) 

But, it is well-known that for the harmonic oscillator (120]) . one has (r 2 ) = (2n + £ + 
3/2)/ (mw). Thus in our case, (r 2 ) = (2n + £ + 3/2) / ^/2~mp~ Q . Using (F54]), we find after 
some algebra that 

(p) « Po- (29) 

This confirms the interpretation of the auxiliary field technique as a mean field 
approximation in which the corresponding field is chosen so that the calculations are as 
easy as possible. 

Using the results developed above, applied to ( jTTl) . one can compute that, for the 
ground state, 



- Sg ' l(A) A + 2 



(30) 



For A = 2, the trivial equality = is found since the method is exact in this case. 
For A = 1, it has been numerically verified that the l.h.s. and r.h.s. of this equation 
are respectively 0.0559 and 0.0525. For A = — 1, these numbers become respectively 
1.25 and 0.764. In both cases, the formula is verified. As expected, the relative error is 
smaller for A = 1 than for A = — 1, since the harmonic potential is "closer" to potential 
r than potential — 1/r. Since E(p ) and E obey the same scaling laws (see section [5]), 
inequality (1301) does not depend on parameters m and a. 
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3.2. Logarithmic potential 

We turn now our attention to the special case of a Schrddinger equation with the 
logarithmic potential {b > 0) 

V(r) = a ln(6r). (31) 

With the choice P(r) = r 2 , equation ([7]) can be easily solved. One finds 




K-\p) 



and (e is the basis of natural logarithms) 

E{p) = x F^-(2n + £ + 3/2) + ^\n 
V m 2 



ab 2 
2~e~p 



(32) 



(33) 



The elimination of p from the minimization condition of this last energy formula finally 
leads to 

ma 2 , , 

Po = 7^—ir-^7^, (34) 



2(2n + £ + 3/2) 2 ' 
a In 



m a 



b(2n + £ + 3/2) 



(35) 



It is well known that the eigenenergies E(m, a, b; n, I) resulting from a Schrddinger 
equation with the potential fl3~Tl) satisfy the property [11] 

E(am, a, b; n, I) = E(m, a, b; n, I) — ^ In a. (36) 

The immediate consequence is that the corresponding spectrum is independent of the 
mass of the particle. It is remarkable that the basic property ( 1361) still holds for our 
approximate expression fl35|) . 

It is worth mentioning that the energy formula ( 135]) can be understood as a 
particular limit case of formula (|25|) . as it was suggested in [12]. The basic idea is 
that the logarithmic potential can be rewritten into a power-law form, i.e. 



In x = lim — (x x — 1). 

A^O A 

Let us then consider the following Hamiltonian 



(37) 



(38) 



The potential term of this Hamiltonian reduces to potential ( !3T|) in the limit A — > 
0. A simple rewriting of formula ( |25l) for a — > a6 A /|A| gives the eigenenergies of 
Hamiltonian fl38l) . We have thus 



E(X) 
and as expected 



2 + A\ (ab x )^ 



2A 



A 

m \+2 



(2n + £ + 3/2)— - 



A 5 



lim E(X) = a In 



m a 



■b{2n + 1 + 3/2) 



(39) 



(40) 
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that is precisely formula (155]) . directly obtained from the logarithmic potential. This 
confirms the idea that the logarithmic potential can be seen as the limit of a power-law 
potential r A when A goes to zero. 



4. The case P{r) 



-1/r 



4-1. Power-law potentials 

Another obvious case of analytically solvable radial potential is the Coulomb potential 



V c (r) = - 



(41) 



whose eigenfunctions and eigenenergies respectively read (see for example [TJ problem 
67]) 

R c (r) = N n ^e^L 2 n ^(2 1 r), E c = ™f + ^ (42) 

with 7 = run/ (n + £ + 1). 

It is thus of interest to choose P(r) = —1/r and to perform the same calculations 
as the ones we did in the previous section. Consequently, we will give less details and 
only mention the main steps of the computations. Equation (jjTj) takes the form 



p = a\X\ r 



A+l 



K-\p) 



P 



a A 



A+l 



(43) 



and the energy spectrum is analytically known once the auxiliary field p is introduced 



as m 



Po 



E C ( Pc 



— {a\\\)*+i (n + 
m 

2 + X\ (o|A|)W5 
~2X 



A 

m A + 2 



1) ; 



n + 



A+l 
A+2 



, 2A 
+ l)A+2. 



(44) 



(45) 



Surprisingly, formula (1431) is identical to (125|) . but the (n + £ + 1) factor due to the use 
of the Coulomb potential P(r) is now present. The (2n + £ + 3/2) or (n + £+ 1) factor is 
thus a "remembrance" of the harmonic or Coulomb approximation. However, the same 
qualitative behaviors in £~ 2 and n~ 2 are obtained. In particular, 



£°(Po)| A=2 




ma 



(46) 
(47) 



- 1 2(n + £+l) 2 ' 

As expected, this method is exact for the Coulomb potential. If, in the formula for 
A = 2, we substitute n + 1 by 2n + 3/2 we recover the exact result. Once more, it is 
remarkable that using a potential P(r) which decreases at large r, one can mimic the 
behavior of the eigenenergies due to a potential V(r) which strongly increases at large 
r. It is worth mentioning that, for large n and £ values, both E H (p ) and E c (p ) are 
proportional to n 2X ^ x+2 ^ and £ 2A /( A + 2 ). Such an asymptotic behavior is in agreement 
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with the WKB analysis which is performed in [TT]. The same behavior is also found in 

The interpretation of the auxiliary field method as a particular mean field 
approximation would be strengthened if we had (p) ~ p as m the previous section, 
with p and po given by ( H3l) and (jSj) respectively. The average value should actually be 
computed with the eigenstates of Hamiltonian (Q in the particular case we are dealing 
with. These eigenstates are actually eigenstates of the following Hamiltonian 

Ha = f - ft, (48) 

2m r 
whose eigenenergies read 

mpl 
2{n + t + Y) 1 

By applying the virial theorem to this last Hamiltonian, one obtains Ea = 
— (po/2) (1/V) ~ — Po/(2(r)). Thus, one has 

<r> « (50, 

and 

(p) = a |A| <r A+1 ) wa|A| (r) A+1 . (51) 

A straightforward calculation shows that these two relations together with the 
definition (jUJ) leads to (p) ~ po, as expected. 

Using the results developed above, applied to (|T7|) . one can compute that, for the 
ground state, 

T(A + 3)" 



E(po) -E 2 

- sgn(A) AT^ 



1 ii • (52) 

2 A+l v 7 

For A = — 1, the trivial equality = is found since the method is exact in this case. 
For A = 1, it has been numerically verified that the l.h.s. and r.h.s. of this equation 
are respectively —0.237 and —0.333. For A = 2, these numbers become respectively 
—1/2 and —1. In both cases, the formula is verified. As expected, the relative error is 
smaller for A = 1 than for A = 2, since the Coulomb potential is "closer" to potential 
r than potential r 2 . Since E(po) and E obey the same scaling laws (see section EJ), 
inequality fl52l) does not depend on parameters m and a. 

4-2. Logarithmic potential 

The resolution of the Schrddinger equation with the logarithmic potential fl3Tl) is even 
simpler with the present choice of P(r). Indeed, one finds in this case 

r = K-\p) = P ) (53) 
a 



and 

F,(n\ = - 

2(n + £+ l) 2 



£(PH- '"f, „, +aln( b -£±). (54) 
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The elimination of p finally leads to 

a(n + £ + l) 2 



m 



a In 



m a 



b{n + £ + l) 



(55) 
(56) 



Again, formulae (135]) and (|56|) exhibit the same formal expression but a (n+£ + l) factor 
is now present. Consequently, the properties mentioned above concerning the spectrum 
and the limit of a power-law potential when A — > still hold in this case. 



5. Scaling laws 



Let us consider the Schrddinger equation 

E(m, 7, a) 



P 



— + 1 V{ar) 



0. 



(57) 



where 7 and a are two real parameters. This last equation can be rewritten with a 
different set of parameters (denoted as m', 7', and a') and with the change of variables 
r = (3 q. After a further multiplication by the constant Xi we obtain 



X: 



Pa 



*(qr) = 0. 



(5f 



'2m'/3 2 

It is readily checked that, with (3 = a /a' and x — /3 2 m'/m, (1581) reduces to (J5j 
provided that the following scaling laws are assumed 

2 m' 
m ' 

(60) 



7 



7 



a 



a' 



(59) 



E(m, 7 ,a)= - —E{m'n',a'). 
\Q :'/ m 



This is the expression of the scaling law for the most general potential V(r). However, in 
the special case of power-law potentials V(r) = ar x , we have a more powerful expression 
due to the fact that we have some freedom on 7 and a parameters. Indeed, if these 
parameters are chosen with the constraint 7« A = a they lead to the same potential. 
This property allows to obtain the absolute dependence of the energies in term of m 
and a parameters. Explicitly, we can write 

E(m, a; A, n, I) = 2^a^m"^ e (A, n, t), (61) 

where e(A,n, £) is an eigenvalue of the dimensionless Schrddinger equation 
a 2 

H = y + sgn(A) I a; |\ (62) 

This Hamiltonian corresponds to the particular values m = 2 and a = 1. This choice 
is made such that, in the exact cases, e(2,n,£) = 2n + £ + 3/2 and e(— l,n,£) = 
— (n + £ + 1)~ 2 . One sees that equations ff25l) and fj45|) resulting from our treatment 
fulfills exactly this scaling law f[6T|) . and that the dimensioned factor gives the correct 
dimension to the energy E(m, a; A, n, £). 
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So an approximate expression for e(A, n, £), that is deduced from the auxiliary field 
technique, can be written as 

(A, n, £) = (*±*) (iv(-PP) („, £)) ^ , (63) 

with the values N^ H \n, £) = 2n + I + 3/2 for the approximation based on the harmonic 
function P(r) = r 2 and N^(n,£) = n + I + 1 for the approximation based on the 
Coulomb function P(r) = —1/r. In the next section we will show how to improve the 
results suggested by the previous studies. 

6. Improved formulae 

Following the form chosen for P(r), two different approximations (163]) are obtained 
for both power-law and logarithmic potentials. A comparison with exact results 
(numerically computed) shows that these formulae are satisfactory but it is possible 
to improve the results. Using scaling laws presented above, it is sufficient to focus our 
attention on the e ( a PP) (A, n, I) quantity. 

We propose here the following expression 

2 

e W(A, n, £) = (^) ^ [6(A) n + £ + c(A)]^ . (64) 
V 2A J 2~ 

Since this formula is exact for A = — 1 and A = 2, we must have £>(— 1) = 1, c(— 1) = 1, 
6(2) = 2, and c(2) = 3/2. Moreover, we can expect that the use of P(r) = —1/r will 
give better results for negative values of A, while the use of P(r) = r 2 will give better 
results for positive values of A. So a natural choice would be 

6(A) = 2, c(A) = 3/2, for A > 0, 

6(A) = 1, c(A) = 1, for - 2 < A < 0. (65) 

The variations of the energy being smooth for variations of A, we can refine our 
assumptions and suppose that 6(A) and c(A) are continuous functions of this parameter 
(the particular case A = will be discussed with the logarithmic potential). So one can 
first think to a linear interpolation between the two analytical cases 

6(A) = *±i c(A) = A±Z. (66) 

Nevertheless, it is possible to greatly improve the accuracy of the approxima- 
tion (164")) by looking at numerical results. Very accurate eigenvalues e n um(A,n, £) of 
Hamiltonian ( |62l) can be obtained numerically with the Lagrange mesh method [14J. In 
order to find the best possible values for the coefficients 6(A) and c(A), we will use the 
following measure 

3 3 

X(A) = ( e num(A,n,£) - e ^(A, n, £)) 2 . (67) 

ra=0 1=0 
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Figure 1. Coefficient 6(A) from (|64p as a function of A, obtained from various 
approximations: numerical fit minimizing the function x{\) (solid line), linear fit 
from |66|) (long-dashed line), hyperbola from (|69|) (short-dashed line), hyperbola 
from ([71]) (dashed-dotted line). This last curve is nearly indistinguishable from the 
solid line. 




Other choices are possible but we find this one very convenient. The coefficients 6(A) 
and c(A) minimizing the quantity xW have been computed for several values of A (see 
figures [1] and [2]) . In both cases, the curve is similar to a section of conic. 
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To determine analytically the best conies fitting the coefficients, it is necessary to 
obtain a third known point on each curve. Solutions of the case A = 1 are known for i = 
in terms of zeros of the Airy function. Approximate eigenvalues of the Hamiltonian ( 1621 ) 
with a linear potential are given by [8j [15] (this is an asymptotic form valid for large n) 

*Ai(l>M)=(^ T J + . (68) 

Equating relations fl68l) and fl64l . we find 6(1) = ir/^/S and c(l) = v^3tt/4. Assuming 
that this procedure can be extended for any values of £, the conies can be determined 
by using the exact values b(— 1) = 1, c(— 1) = 1, 6(2) = 2, and c(2) = 3/2. We find 
that hyperbolae give the best results (minimal value of xW)- The coefficients are then 
given by 

a = V37T, 

(4a - 18)A + (18 - 2a) 



6(A) 



(3a- 15) A + (21 -3a)' 



(7a-36)A+(36-5a) 
C(A) = (6a-32)A+(40-6a)- (69) 
These relations could look quite complicated but it is worth noting the symmetries 
existing between the absolute values of the numbers present in the formula. When 
a number appears twice in a coefficient of the numerator (denominator), it is the 
arithmetic mean of two other corresponding numbers in the denominator (numerator). 
A number which appears twice in c(A) is the double that the corresponding number in 
6(A). The quality of this approximation can be estimated in figures [1] and [2], and in 
tabled] where results obtained from the following formulae (computed in the framework 
of the WKB method) are also given [TT] 

CWKB(W)= + f Av ^r(|) + ^ ( 2n + £ + ^) 1 for A>0 ' 




lAlVTFr(i-i) 

for - 2 < A < 0. (70) 
It is clear that the approximations are mainly good for —1 < A < 2. Let us note that 
e WKB(l' n ' 0) = e Ai(l' 0)- Thus, approximations (TTD]) and (IMj) with (|69l) give the 
same results in the cases A = 2, —1 (exact results), and A = 1 for I = (approximations 
that tend to the exact results for large n); however they differ for all other sets (A, n, i). 
Our formulation has the advantage to deal with a unique consistent formula, valid for 
any A, whereas the WKB approximation requires two different expressions, one for A > 
and one for A < 0. A comparison with exact results (see tableland table [2]) shows that 
our formula is better than the WKB formula in any case. 

Other approximate forms for the eigenvalues of Hamiltonian ( 1621) exist and give 
good results [13l [16], but they are not considered here because they do not give the 
exact result for the cases A = — 1 and A = 2. 
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Table 1. Values of xW defined by (|67|) . as a function of A, for various approximations. 
x{— 1) = x(2) = in all cases. The value A = corresponds to the logarithmic potential 
(see text). 
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It is possible to improve the value of the quantity x(A) by fitting the hyperbolae 6(A) 
and c(A) directly with the numerical values enum(A,n, £). Several choices are possibles. 
We give here a parameterization with simple numbers (below 100) which give good 
results 

L/x . 41A + 86 /x . 5A + 17 

6(A) = 13AW C(A) = 2ATI4- (71) 
The quality of this fit can be appraised in figures [1] and El and in table [Q Obviously 
the numbers in ( |7T|) depend on the points used for the fit but also on the particular 
choice of the function xW- Other definitions, relative error instead of absolute 
error or different summations on quantum numbers, would have given other numbers. 
Parameterization (ITTT) gives better global results than parameterization fl69|) . except for 
small values of £ and large values of n. This behavior is expected since formulae ( |69l ) 
rely on the formula (168]) whose accuracy increases with n |15j . 

In order to have a better estimate of the qualities of our results, we compare in 
table El the dimensionless energies e(A = l,n,£) obtained from various approximations. 
One can remark that equation (IS~4"|) with approximation (|69|) gives good results, as good 
as or better than e\\TKB' anc ^ that ec L ua ti° n (E1D with approximation ( ITTi) gives very 
good overall results, much better than the previous ones for i ^ 0. This is particularly 
significant for large values of n and I quantum numbers. These conclusions remain valid 
for any value of A. 

Let us now consider the dimensionless logarithmic potential 
q 2 

H =— + \n\x\. (72) 
4 1 1 
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Table 2. Values of e(A = l,n,£) for various approximations. For each set (n,£), 
first line is the result obtained by numerical integration [enum(b n, £)], second line 
by equations (|64|) and (|69|) . third line by equations (j64|) and (|7T|) . and fourth line by 
equation g£jj) [e WKB (l, n, I)]. 



t e(l,0,l) e(l,l,l) e(l,2,l) e(l, 3, £) 

1.47292 2.57525 3.47773 4.27536 
1.46167 2.57138 3.47560 4.27394 
1.47214 2.56575 3.45909 4.24854 
1.46167 2.57138 3.47560 4.27394 

1 2.11746 3.07701 3.91056 4.66528 
2.11057 3.08645 3.92585 4.68320 
2.11929 3.08132 3.91032 4.65894 
2.05470 3.04039 3.88505 4.64587 

2 2.67619 3.54649 4.32712 5.04580 
2.67098 3.56156 4.35159 5.07529 
2.67874 3.55678 4.33684 5.05198 
2.57138 3.47560 4.27394 5.00345 

3 3.18188 3.98898 4.72763 5.41584 
3.17757 4.00682 4.75742 5.45277 
3.18469 4.00231 4.74332 5.43029 
3.04039 3.88505 4.64587 5.34868 



Using the results above as well as (HOI , it is clear that a good approximate form for the 
eigenvalues is 



The quality of this formula can be estimated in table [U 
7. Conclusions 

In this paper, we propose a new method to obtain approximate expressions for the 
eigenenergies of a Schrddinger equation based on any local central potential V(r) and 
valid for any radial n and orbital I quantum numbers. The method relies on the 
introduction of an auxiliary field p. Given a potential Pir) for which the analytical 
spectrum of the Schrddinger equation is known, we presented an algorithm which 
transforms the solvable Hamiltonian into the Hamiltonian which contains the original 
potential V(r) and for which the eigenenergies are not known analytically. This 
transformation needs an auxiliary field which makes a bridge between V(r) and P(r). 
If the auxiliary field is considered as an operator p which cancels the variation of the 
transformed Hamiltonian, the algorithm gives the exact solution. We propose to start 
from the solvable Hamiltonian and the corresponding eigenenergies expressed in terms 




(73) 
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of the auxiliary field, and consider this field not as an operator but as a pure number. If 
this number is chosen by a minimization procedure, the resulting energies gives a good 
approximation to the exact values. If P(r) = V(r) our method gives the exact results. 

This method is applied to a power-law type of potential with power A. The resulting 
expression for the energies satisfies the scaling law prescriptions and gives the good 
asymptotic behavior for the radial and orbital quantum numbers. Comparing the 
results coming from the use of P(r) = r 2 and P(r) = — 1/r, one sees that they are 
identical provided that the (2n + / + 3/2) term, remembrance of the harmonic oscillator, 
is replaced by (n + I + 1), remembrance of the Coulomb potential. Owing to this 
remarkable property, we proposed new analytical formulae, valid for any A, n, £, that 
give results much better than the equivalent ones which could be found in literature 
since they are able to reproduce the exact results with an accuracy often better than 
10 -3 . These formulae can be used for example if one wishes to determine the parameters 
of a power-law potential that fits the experimental data as better as possible. We also 
showed that the results for logarithmic potentials are obtained very easily from those of 
the power-law potential if one takes the limit A — ► 0. 

The main interest of our method is to obtain analytical formulae. Its main difficulty 
consists of the inversion of some algebraical equations. In the WKB method, analytical 
formulae can also be obtained provided the integral over the radial momentum can be 
performed. Depending on the particular problem considered, one or the other method 
can be more efficient. So the new method based on an auxiliary field, proposed here, 
can be a useful supplementary technique to solve eigenvalue problems. It seems very 
promising to obtain analytical expressions for the eigenenergies of a Schrddinger equation 
with a number of interesting potentials. Some of them are under study. Analytical 
approximate eigenvectors are also given by this method. Their study will be considered 
in a subsequent paper. 
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